Contribution of One-Time Pair Correlation Function 
to Kinetic Phenomena in Nonequilibrium Gas 

A. Ya. Shul'man 
Institute of Radioengineering and Electronics of the RAS 
Moscow 101999, Russia. E-mail: asli@cplire.ru 

Abstract 

It has been established in nineteen seventies that in nonequilibrium case the pair colhsions generate 
non-zero two-particle correlations which are non-diagonal in momentum space and give the essential 
contribution to the current fluctuations of hot electrons. It is shown here that this correlations 
give also a contribution to the collision integral, i.e., to kinetic properties of nonequilibrium gas. 
The expression for electron energy loss rate P via phonons is re-derived in detail from this point of 
view. The contribution of the non-diagonal part of the nonequilibrium pair correlator to phonon- 
electron collision integral and to P is obtained and explicitly calculated in the electron temperature 
approximation. It is shown that these results can be obtained from stochastic non-linear kinetic 
equation with Langevin fluctuation force. Such an approach allows to formulate the simple general 
conditions under that a contribution of two-particle correlations might be essential in kinetics. The 
contribution obtained does not contain the extra powers of small gas parameter unlike the equilibrium 
virial decompositions. 

1 Introduction 

The derivation of kinetic equation for one-particle distribution function F{p, r, t) from BBGKY chain 
leads to the collision integral expressed in terms of the two-particle correlation function. The latter is in 
turn expressed by F{p,r,t) using Boltzmann's StoBzahlansatz that implies the neglect of initial one-time 
correlation function g2 {pi,p2,r,t) , which is non-diagonal in the momentum space. In nonequilibrium 
case, however, it has been established 0] — that pair collisions generate the non-diagonal two-particle 
correlator that gives the essential contribution to the current fluctuations of hot electrons. Here we 
show that this correlator gives also a contribution to the collision integral, i.e., to kinetic properties of 
nonequilibrium electrons. The expression for electron energy loss rate P via phonons is re-derived in 
detail from this point of view. 

In fact, the formula for P in terms of electron density-density correlation function {pkpk)^ was 
obtained by Kogan using Fermi Golden Rule. However, to come to the well-known result for P in 
the electron temperature approximation the fluctuation-dissipation theorem was implicated to express 
(pkPk)^ in terms of the dielectric function of the electron gas. This step is defensible in the case 
of the true equilibrium electron gas only and it gives rise to results equivalent to the StoBzahlansatz. 
In this report the expression for P in terms of {pkpt)^^ is derived from quantum kinetic equation for 
phonons interacting with nonequilibrium electron gas. The contribution of the non-diagonal part of the 
nonequilibrium pair correlator to phonon-electron collision integral and to P is obtained and calculated 
explicitly in the electron temperature approximation. 

It is shown that these results may be obtained from stochastic non-linear kinetic equation with 
Langevin fluctuation force. Such an approach allows to formulate the simple general conditions when 
the two-particle correlator contribution might be essential in kinetics. The derivation of this equation 
based on the procedure of partial averaging over microscopic-scale fluctuations similar to the procedure 
usmg by H. Mori [| to derive the generalized classical Boltzmann equation for dense gases and liquids. 

'^Based on a poster presented at the conference "Progress in Nonequilibrium Green's Functions, Dresden, Germany, 
19.-22. August 2002" 
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However, in our case the contribution obtained does not contain the extra powers of small gas parameter 
of the kinetic theory. 



2 Kinetic equation for phonons 

2.1 Basic designations 

Hamiltonian H = Hp + Hph + Veph 

Electron part = ^ hilka^Uk + Vee where Vee stands for the electron-electron interaction energy 

k 

Phonon part Hpfi — ^ hujfbtbf 
f 

Electron-phonon interaction 

Veph = ^ E (c/P/ fe/ + c}Pfb^) 

Spacial Fourier components of the electron density operator 

k ^ 

Phonon distribution function 

m) = {(bpf)) 

Electron distribution function 

The angle brackets denote the full averaging with the density matrix corresponding to the nonin- 
teracting electrons and phonons at an initial instant of time. The subscript t indicates the Heisenberg 
representation for electron and phonon operators. 



2.2 Phonon-electron collisional integral 

and the electron density-density correlation function 

We obtain the kinetic equation for slow variables Nf{t), smoothing Heisenberg equations of motion over 
microscopically-large time interval At. Such a technique was suggested by N.N. Bogolyubov and N.M. 
Krylov (see the book ||], pp. 5-76) and it has been used in |^ for the derivation of the linearized kinetic 
equation for the occupation number fluctuations with Langevin sources of fluctuations. With accuracy 
to the second-order terms in e ~ ph interaction we have 

, t+At u , 

- [s+{t) J dt'J dt" [V^p,{t'% [Vp%it'), (b+bf),]] S{t)j (1) 

The superscript denotes operators in the interaction representation with respect to the electron-phonon 
interaction Vep- 

Let us remind interrelations between S'-matrix and the time evolution of any operator O in the 
Heisenberg and interaction representations: 

Uit,to) = e-*^(*-'«) = U"{t,to)S{t,to), U%t,to) = e-*(^'=+^'"')(*-*°\ 

d{t) = U+it,to)dito)U{t,to) = 5+(t,to)0"(i,to)^(i,to) (2) 
ihdS/dt^V,l^,{t)S, S{to) = 1, O°(t,to) = t/°+(t,io)0(to)t/°(i,to) 

Calculations of the internal commutators in Eq.([l]) lead to the desired result involving the electron 
density-density correlation function: 



where the commutativity of free electron and phonon operators in Eq.(|l|) and the condition uifAt ^ 1 
have been used. 
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The function 

oo 

(p/P/ )^ = y" dr (^pfit + T)p+{t)) exp {iuT) (4) 

— oo 

is the spectral density of the correlation function of the electron density fluctuations. It is interested to 
remark that the right hand side of Eq. (0) may be rather interpreted as the emission and the absorption 

of phonons by the electron fluid with the structure factor \PfP^j than as individual electron- phonon 

collisions. 

2.3 Spectral density of the electron fluctuations 

The sought spectral density must be calculated without Vep taken into account. From equations of motion 
for the electron operators we obtain in the random phase approximation: 

|e(w,/)| 

where e((jj, /) is the dielectric function of the electron gas and 
PfP^T^ = E {Pf{k)p^{ki) 

k,ki 



is the spectral density of the density-density correlation function for non-interacting electrons. The 
one-time correlation function is 

{pfik)p+{h))^ EE a,+ /a+ / )^ . (6) 

The one-time correlation function of the non-interacting electrons is independent of time (in steady state) 
under condition that the electron energies hilk — fi^k'^ /2m obey the equality 



As long as / 7^ it is necessary to he ki = k . Therefore, this condition provides the one-time 
correlator to be the integral of motion during the time interval At like the distribution function Fk . 

2.4 Kinetic equation with account for the contribution of nonequilibrium 
pair correlations 

The substitution of the explicit expression for correlators to the kinetic equation (H) gives rise to 



dN 



L 

dt 



voh?\e{uiiJ)V fe-i fc+i fe+i ^-i/t \ ''+^ ''/ 
Taking into account the commutation rules a^^ + ak^ ct^^ = 1 it is easy to obtain 



3 



According to the results from Gantsevich et al. ||T| 

{ata+^ak,ak) ^ FkFk,+<j){k,ki) (10) 

where (jj (k,ki) is a noncquilibrium correction to the pair correlator. 

Finally, the kinetic equation for phonons with corrections for the non-equilibrium part of the electron 
pair correlator takes the form: 

dNf ^ 
dt 

Z^lSjr 5 (1 - - ^ - ^ - ^'-f - - (11) 

It should be stressed that for ordinary classical equilibrium electron gas there are no pair correlations 
between different states in momentum space, i.e., 4> (k, fci) = 0. This fact can be easily seen from the form 
of the Gibbs distribution with the usual kinetic energy part of the Hamilton's function K — J^i P?/2to 
since there are no terms proportional to PiPj stand for the particle enumeration). 

3 Interrelation between non-linear stochastic kinetic equation 
and correlation corrections 

Following Q — and using partial averaging over fast and short-wave fluctuations under conditions that 
the slow variables like Nf{t) and Fk{t) are constant instead of the full averaging over all states of the 
system under consideration, we find the non-linear stochastic kinetic equation with Langevin sources of 
fluctuations in the form 



X 



{[Nf + [Pk (l - Fk-f)] S {ojk - ujk-f - c.) - (12) 

-^/E [Pk {l^ Fk+jy\5{ujk-0Jk+s 

The similar stochastic kinetic equation can be derived for the fluctuating electron occupation numbers 
Fk with accounting for Vee and Veph interactions . Both equations are nonlinear in and bilinear in 
NfFk- Here tilde denotes that the slow fluctuations are still kept. Langevin source of fluctuations can 
be expressed as the difference between the exact time derivative of the phonon occupation number and 
its approximation 

t+At p 1 t+At p 1 

Jf{t) = j^^S^{t) J dt' K%it'),N^fit') S{t)~j^^S+it) J dt' V^^,{t'),NOit') Sit) 

t -' t '- -' 

It is seen that the partial average of Langevin sources themselves is equal to zero Jf{t) — . The 

partial correlation function i[Jy(ii), Jq{t2)]+ is expressed by means of the transitions' rates involved in 
the collision integral similarly to the case of the linear stochastic kinetic equation for fluctuations of the 
occupation numbers 

Taking full average from stochastic kinetic equation and using the evident properties (see Eq.([lO|)) 



Fkj = Fk; (^FkFk.j - FkFk,+^{k,ki),ky^ ki,and also (J/(t)) = (13) 

we can easy obtain the same kinetic equation for phonon occupation numbers with correlation corrections 
that was derived above. Moreover, from this point of view we can immediately write the respective terms 
in the Boltzmann kinetic equation for neutral particles and bring up the question about a possible 
contribution of nonequilibrium particle-particle correlations to the equations describing phenomena like 
flow turbulence. The similar statement concerns to electron-phonon correlations described by expressions 
like (SFkSNf ) in the case when a phonon nonequilibrium takes place. 



f 

In conclusion of these general considerations it needs to underline that the contribution of non- 
equilibrium correlations to kinetics should be expected in that cases when the excitations considered are 
scattered by nonequilibrium many-particle system. 
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4 Electron energy loss rate via phonons 



The energy loss rate of electrons due to phonons per unit volume can be obtained using the relation: 



1 dNf 

-E^/^- (14) 



vo ^ at 

It is seen from the kinetic equation for phonons Eq.(^l]) that the loss rate can be broken up into two 
parts: 

P = Pl+P2, (15) 

where on the right hand side the first term is related to the one-particle electron distribution function 
and it is described by the known expressions. The second one is due to the correlation contributions 
determined by the function </>. We assume during the calculations of P that the phonons are in thermal 
equilibrium and the electron-electron collisions are sufficiently frequent in order to provide the use of the 
electron temperature approximation for hot electrons. As a result of these assumptions the following 
formulae may be derived: 

P,{T) = y ^HiOf [N%T) - NA 



X 

k 



and 



2 

P2(T, To)^-y fj^f hiof J2^ik + f,k)S (n^+f - Qk - iof) . (17) 

Here Nj{T) = [exp (fiwy/T) — 1] ^ is the equilibrium Bose-Einstein distribution with hot-electron 
temperature T and Nf is the true noncquilibrium phonon distribution function. Further, we will consider 
the equilibrium phonons, i.e., Nf — N*^- (Tq). In this case the expression (|l6|) for Pi coincides with well- 
known Kogan's formula except for dielectric function in the denominator describing the dynamic screening 
of the electron-phonon interaction. 

The function </> (fc, fci) was found in |Q under accepted assumptions and can be presented in the form: 

cj>ik,h) = Q{T,To)^^, (18) 
where F^{T) is the Maxwell distribution of non-degenerate hot electrons and 



rp2 

Q{T,To) = 

nc„ 



p(r,To)/(r-To) 

dP/dT 
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(19) 



The expression for in g was obtained without accounting for the contribution of </> itself in P . 
Therefore, the quantity Pi was only implied as the total energy loss rate P. Now evidently we have the 
equation for determination of P2(T, Tq) since this quantity is defined by the right hand side of the equality 
(^) where <j){k + /, k) depends itself on P2(r, To) accordingly to Eqs.(|l|),(|5Hl|). The situation is very 
resembling the self-consistent field approach in the plasma theory. Now we have to solve simultaneously 
the set of kinetic equations for the one-particle distribution function and the pair correlator. 

It needs to stress that the substitution of the total loss rate P instead of Pi in the expression (|l^) for 
(f> {k, ki) is justified because all symmetry properties of the electron-electron and electron-phonon integrals 
of collisions used in during the derivation are not destroyed by introducing (j){k,ki) in collision integrals 
for electrons and phonons. 



It is of importance to note in conclusion that the expression (18) for (j){k^ki) found in the limit 
of high frequency of the electron-electron collisions does contain coupling constants of neither electron- 
electron nor electron-phonon interactions. Therefore, the corresponding corrections to the integral of 
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phonon-electron collisions do not contain additional powers of the small parameters of the kinetic theory. 
At low heating it is proportional to T — Tq like Pi{T,Tq), as it can be obtained from the expression 
(|l9| ) for Q(r, To). However, the solution ( |l8|) for (f>{k,ki) has been obtained in ||] in the case of the 
non-degenerate electron gas. Hence Fj}{T) <C 1 and the product of two distribution functions in Eq.dl^) 
makes P2 <C Pi in the case of scattering electrons by acoustic phonons since hcoj <C T, Tq and Nj ^ 1 
for phonons involved. The situation can change in the case of low-temperature scattering of electrons by 
optical phonons under condition fkoo ^ T,Tq. The 4>{k,ki) corrections should be more pronounced for 
the hot-electron phenomena in the degenerate electron gas of semiconductors or metals [^. 
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